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Abstract. The large N limit of mean spectral density for the ensemble of N x N 
sub-unitary matrices derived by Wei and Fyodorov (J. Phys. A: Math. Theor. 41 
(2008) 50201, Ref. p]) is calculated by a modification of the saddle point method. It 
is shown that the result coincides with the one obtained within the free probability 
theory by Haagerup and Larsen (J. Funct. Anal. 176 (2000) 331, Ref. [5]). 
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Introduction 

In many different physical and mathematical problems (see e.g. p] and references 
therein) it is of interest to know the spectral density of sub-unitary ensemble of random 
matrices of the form 

T = UH (1) 

where U is N x N unitary matrix distributed according to the standard Haar measure 
and if is a diagonal matrix of the form 

Hij = S i:j y/gl (2) 

where each gi obeys < Qi < 1. 

The exact formula for the spectral density of these matrices for finite N and 
arbitrary gi has been obtained in pQ by a supersymmetry method. In the conclusion 
of that paper the authors wrote: "An interesting problem would be to investigate the 
density of complex eigenvalues in the limit N — > oo assuming that g has a finite density 
v(g) = jj J2i $(.9 ~~ 9i) of eigenvalues in an interval of the g-axis". 

First of all we would like to stress that such type of asymptotic questions is very 
natural within the free probability theory [2]-@]- In fact, an explicit solution to this 
problem has been found in [5]. The purpose of this note is to demonstrate how this 
answer can be derived from the exact formulas obtained in pp. 
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Free probability solution 

For general the matrix (OQ) have complex eigenvalues X a . The question we address is 
to calculate the mean density of these eigenvalues in the limit of large matrix dimensions 

v T {z)= lim hj^Siz-X,)). (3) 

i 

Here (...) denoted the average over the Haar measure of matrices U. 
Assume that limiting positive moments of gi in (T5]) 

1 N 

i=l 

are finite which is equivalent to the existence of finite limiting measure for g^. 

According to the theorem 4.4 of [5] the limiting measure vt{z) ([3]) of ensemble of 
matrices (CQ) is non-zero in the annulus 

(^(-1)) _1 < N 2 <M S (1)- (5) 

If Hg(—1) diverges the lower bound is set to zero. 

Inside this annulus the mean number of eigenvalues of matrix flTJ in complex plane 
in polar coordinates z = re 1 ^ is 

dn = u T (r)rdrd(p , (6) 

and the integrated spectral density is determined from the equation 
rHv) 

y = UT(r)rdr . (7) 

Jo 

The function F(y) here is calculated as follows. Define first the generation function of 
integer moments ^ g (u) by the formula 

1 N 

VJu) = lim - V^- (8) 
Then one has to invert this function, i.e. to find the function x = Xg(y) such that 

v g (x(y)) = y- (9) 

The knowledge of this function permits to calculate the so-called ^-transform of g which 
is typical in free probability theory [3] 

s g (w) = ^^x a M . (io) 

Finally the function F(y) which enters to the mean density of eigenvalues (jTJ) takes the 
form 

F(y) = (S(y-l)y 1/2 . (11) 
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In other words, the integrated mean density of eigenvalues, y = y(r), is the inverse 
function to F(y) i.e. r = F(y). From the above equations one easily gets that y(r) has 
to be determined from the equation 



yr i 



y-i- 



(12) 



For clarity let us consider the simplest example of truncated random unitary matrix 
ensemble investigated in [6]. It corresponds to the case when g% = \ for i = 1, . . . , M 
with all others gi — and M/N — > \i when N — > oo. In this case a = fi and 6 = 0. The 
function \l/ in (jSj) is now 

u 



0: 



(13) 



1 - u 

from which it follows that the solution of (fl2l) is 
1 — /i 

f = 1 2 

1 — 

which corresponds to the radial spectral density equal to 

1 - /i 
1 ' tt(1 -r 2 ) 2 

This result agrees with Eq. (19) of [6] up to the normalisation. In that paper the authors 
normalized the density to the number of non-zero eigenvalues, but in the above formulas 
the density is normalized to the total number of eigenvalues which gives an additional 
factor of ii. 



(14) 



(15) 



Derivation of asymptotic formula from exact solution 

Assume that gi are ordered: < g% < g% . . . < g^- Then according to Eq. (2.3) of [1] 
the exact mean density can be written as follows 







\z\ < 9i 



v T (z 



jj Eifc+i F A (gi) gk < \z\ 2 < g k +i ■ 




9n < \z\ 



Here Fjx(gi) is the following function 



F A {gi 



(gi 



\2\N-2 /"oo 



Ujd^-dj) Jo (1 + *) 



Ndt 



\N+2 



n 



N-t+ f^(Nt- 1) 



The density in ffl6l) is normalized in such a way that 



v T (z)d\z\ 



1 . 



For the later use it is convenient to change t — > u as t 



u)/u. It gives 



N 



F A {gi 



w(u)du 



N 



(l-u)/u + gi /\z\ 2 
l + (l-u)g t /(u\z\i) 



(16) 



(17) 



(19) 
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Figure 1. Contour of integration in <P5|) . v p = g p /(g p —\z\ 2 ) are poles of the integrand 
in (p?5|) . The positions of the fixed saddle point v — y obeying ([25)1 and of the moving 
pole v = u are indicated by the arrows. 



where 



W (U) 



u 



N 



N 

n 



u 



u\z 



1 



(20) 



As this function is symmetric functions of g^, the sum entered to ffT6l) can be represented 
as the contour integral. Define the function 

{q- \z\ 2 ) N ~ 2 



f(q) 



Uf=i(q-9j) 



It is clear that 



u T {z 



Z7T1 



N 



u)/u + q/\z\ 



(21) 



(22) 



1 + (1 -u)q/(u\z\ 2 )_ 

where the contour Ck enclosed poles gk+i, gk+2, ■ ■ ■ ,9n- 

Changing the variable q = v \z\ 2 /(v — 1), one gets the convenient representation of 
the exact eigenvalue density (|T6l) 

1 f dv 



v T \z\ 



27ri M 2 fa w(v) J 



w{u)du 



N 



u + v — 1 



(23) 



and the contour C' k encircles the images of the same poles as above but has to avoid the 
pole at v = u with < u < 1 as indicated at figure [Q Notice that the functions in the 
denominator and in the numerator are the same which reveals the supersymmetry-like 
structure of these integrals. 

The function w(v) in ( 1201) can be rewritten in the form w(v) = exp[N<&(v)] where 



1 N ( 
$( w )=ln( u )+ JJlnh 

i=i ^ 



v - 1 

v\zr 



(24) 
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At large N the integral over v (as well the integral over u) can be calculated by the 
saddle point method. The saddle point v = y{\z\ 2 ) has to be determined from the 
condition 

N 



i - 1 E- 



(y -i)9j/(y\z\ 



1 



i-y 



y-i 

y\z\ 2 







(25) 



where we introduce the same function ^ g (u 



as in ([8]). From this expression it follows 
that the saddle point obeys the same equation as (fi"2l which justify the notation y for 
the saddle point. 

As the same saddle point appears also in the integration over u, it will be important 
only when it lies between the limits of integration over u, i.e. when it obeys the 
inequalities < y < 1. It is straightforward to conclude that it will be the case 
only when the inequalities ([5]) for \z\ 2 are fulfilled. From these arguments it follows that 
when N — > oo the dominant contribution to the integrals (1231) comes from a vicinity of 
the above saddle point y = y(\z\ 2 ). The difference between the integrals over u and over 
v is that in the former the integration is performed along the real axis but in the later 
the integration contour should be deformed to pass through the saddle point parallel to 
the imaginary axis. As the integrand over v decreases quickly at large v the necessary 
deformation of the integration contour is legitimate but as there exists a pole at v — u 
such deformation will result at an additional contribution from this pole (cf. figured]). 
From ( l23i) it is straightforward to check that this contribution takes the form 



(2u - l)du 



1 



■y(y-K 



(26) 



After deforming the contour, the integration in the saddle point approximation is 
standard and the final answer is 



v T (z\ 



N- 



— { wm +yiy ~ 1) 



(27) 



Here $"(?/) is the value of the second derivative of the function ( 1241) calculated at the 
saddle point y obeying (fl2l . From (|25|) one gets 



*"(!/) 



y(y 



l - 



y \z 



2 9 



y 



y\ z \ 



(28) 



where prime indicates the differentiation over argument of the function. One can check 
that $"(?/) is negative as it should be for the application of the saddle point method. 
From (TI2D it follows that 



dy 

d|zp 



1 - 



y \z 



2 9 



y-i 

y\z\ 2 



y 



y\z 



4 9 



y 



y\ z \ 



(29) 



Expressing 



y\z\ 2 



dy 

d\z\' 



through $"(?/) from ff28l) and substituting to ff29l) one gets 



+ y(y 



(30) 
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Taking into account that $"(?/) < one finally concludes that the combination in (1271) 
equals ck//d|;z| 2 , or 

W) — ► ^ (31) 
which exactly corresponds to the free probability result (TT2T) derived in [5]. 
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